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ABSTRACT 

We conducted VLA observations at 0."06 resolution of the 22 GHz water 
masers toward the Class source S106 FIR (d = 600 pc; 15" west of S106-IRS4) 
on two epochs separated by ~3 months. Two compact clusters of the maser 
spots were found in the center of the submillimeter core of S106 FIR. The 
separation of the clusters was ~ 80 mas (48 AU) along PA. = 70° and the size 
of each cluster was ~20 masxlO mas. The western cluster, which had three 
maser components, was 8.0 km s _1 blueshifted and the eastern cluster, which 
has a single component, was 7.0 km s _1 redshifted with respect to the ambient 
cloud velocity. Each component was composed of a few spatially localized maser 
spots and was aligned on a line connecting the clusters. We found relative 
proper motions of the components with ps 30 mas year -1 (18 AU year -1 ) along 
the line. In addition, a series of single-dish observations show that the maser 
components drifted with a radial accelerations of ~ 1 km s -1 year -1 . 

These facts indicate that the masers could be excited by a 10 AU-scale 
jet-like accelerating flow ejected from an assumed protostar located between 
the two clusters. The outflow size traced by the masers was 50 AUx5 AU 
after correction for an inclination angle of 10° which was derived from the 
relative proper motions and radial velocities of the maser components. The 
three-dimensional outflow velocity ranged from 40 to 70 km s -1 assuming 
symmetric proper motions for the blue and red components. Since no distinct 
CO molecular outflows have been detected so far, we suggest that SI 06 FIR 
is an extremely young protostar observed just after the onset of outflowing 
activity. 

Subject headings: ISM: jets and outflows -- ISM: individual (S106 FIR) - 
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1. INTRODUCTION 



The formation mechanism of jets and molecular outflows from young stellar objects 
(YSOs) and their roles in the star formation process are not well understood observationally 
or theoretically. One of the best observational approaches is to investigate very young 
protostars which are in an evolutional stage just after the onset of the outflow activity. 
Andre, Ward-Thompson and Barsony (1993) identified such protostars with sources which 
have not been detected in the near-infrared and which display a blackbody-like spectral 
energy distribution (SED) which peaks in the submillimeter. They proposed to classify 
these objects as 'Class 0' sources, being in a younger evolutional stage than that of 'Class 
I' sources defined by Lada (1987). The cold SEDs of Class sources suggest that the 
bulk of final stellar mass has not yet been accumulated. Class sources are thought to 
be in their main accretion phase and they are known to have very powerful "jet-like" CO 
bipolar molecular outflows which are usually observed with typical sizes of ~ 0.05 — O.lpc 
(e.g., Bachiller et al. 1991; Bachiller et al. 1996; Bachiller 1996; Yu and Chernin 1997). 
In contrast, the CO outflows from Class I sources are poorly collimated and much less 
powerful (Bontemps et al. 1996). It is believed that the acceleration and collimation of 
the outflows involve magnetohydrodynamic processes in the vicinity of the central stars 
(r <^ 1 AU). For example, the outflows are thought to be centrifugally accelerated from 
magnetized Keplerian disks (e.g., Uchida and Shibata 1985; Kudoh and Shibata 1997; 
Pudritz et al. 1991) or thought to be accelerated from the interaction regions between the 
stellar magnetospheres and the surrounding disks (e.g., Shu et al. 1988; Najita et al. 1994). 
Millimeter-interferometer observations with small arrays have extensively investigated the 
kinematics of gas around Class sources. Their spatial resolutions, however, were at most 
arcsecond scale (corresponding to a 100 AU scale in nearby star forming regions). It is 
essential to attain higher spatial resolution - down to subarcsecond-scales - in order to 
understand outflows in the vicinity of Class sources. 
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Interfero metric imaging of the H 2 maser line at 22 GHz with the VLA and VLB A 
provides an excellent tool for studying subarcsecond structure in protostellar jets very close 
to protostars. It has been suggested that the maser lines in low-mass star forming regions 
are associated with outflows, because the luminosities of the masers correlate with mass-loss 
rates derived from CO molecular observations (Felli et al. 1992) and with the luminosities 
of 6 cm free-free continuum emission from ionized, collimated outflows (Wilking et al. 
1994; Meehan et al. 1998). Clear evidence for that was obtained by high spatial resolution 
studies with the VLA-A array (e.g., Chernin 1995) and with the VLBA (e.g., Claussen et 
al. 1998). These observations revealed that the masers originate behind the shocked gas 
very close to the central protostars. Interferometric observations of the maser lines have 
three advantages. First, we can attain high angular resolutions ( ~1-100 mas ) and high 
relative positional accuracies (typically ~ 10% of angular resolutions). Second, the maser 
lines are free from extinction even in deeply embedded sources. Third, we can attain high 
velocity resolution (<^ 0.1 km s^ 1 ) with radio spectroscopy. However, such studies have just 
begun toward the Class 0/1 sources. 

We have found that Class sources provide the best targets for studying gas kinematics 
very close to protostars through observations of the H2O maser lines. Masers have been 
detected for 18 of 30 Class sources (detection rate is ~ 60%) in our ongoing multi-epoch 
survey toward Class 0/1 sources with the Nobeyama 45m-telescope (Furuya et al. 1998). 
On the other hand, we measured a detection rate of only rs 10% for Class I sources. This 
estimate is based upon data compiled from our survey and previous ones (Wouterloot and 
Walmsley 1986; Wilking and Claussen 1987; Cessaroni et al. 1988; Comoretto et al. 1990; 
Terebey et al. 1992; Felli et al. 1992; Palagi et al. 1993; Xiang and Turner 1995; Claussen et 
al. 1996), although sensitivities and angular and velocity resolutions differ. We interpreted 
this result as follows. Since Class sources possess a large amount of circumstellar gas, 
the gas shocked by outflow can easily reach the physical conditions necessary to excite the 
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maser emission (10 6 cm~ 3 <^ n(H 2 ) <^ 10 9 cm~ 3 with T ^ a few hundred K ; Hollenbach 
1997, or n(H 2 ) £ 10 9 cm~ 3 with T ^ 300 K ; Elitzur et al. 1989 ). In addition to providing 
appropriate conditions near the outflow origin, edge-on disks could also enhance the 
detectability of maser emission by maximizing the path length of velocity coherent gas ( 
Elmegreen and Morris 1979). In this case, we could investigate kinematics of protostellar 
disks using the maser probe as shown by Fiebig et al. (1996) in IRAS 00338+6312. 

In this paper, we have studied the 22 GHz H 2 maser emission toward the Class 
source S106 FIR. A compact, bipolar HII region lies near S106 FIR. A distance to S106 
molecular cloud is 600 pc (Eiroa, Elsasser and Lahulla 1979; Staude et al. 1982). The 
exciting star S106-IRS4 (other names are S106-IRS3, S106-IR and S106 PS) whose spectral 
type is O8V-B0V (Gehrz et al. 1982), is located 15" east of S106 FIR. S106 FIR is a Class 
source isolated from this bipolar HE region. A dense core with a diameter of 9000 AU was 
observed around S106 FIR in 450, 800 and 1100/im continuum emission using the JCMT 
(Richer et al. 1993), while no emission at 20 /im was detected (Gehrz et al. 1982). The 
SED of S106 FIR suggests that this source is a Class source. Analysis of the SED implies 
that the luminosity of S106 FIR is in the range of m 24 to 10 3 Lq (Richer et al. 1993). It 
should be noted that a distinct CO molecular outflow has not been detected in S106 FIR, 
making it unique among all Class sources searched (Hayashi et al. 1993; Bachiller 1996; 
Furuya et al. 1999), although there is a possibility that wing emission from the S106FIR 
outflow may be hidden by broad line emission associated with the HII region. The peak of 
the submillimeter core (see Fig. 4 of Richer et al. 1993) corresponds to the position of the 
H 2 masers. The maser emission was first detected by Stutzki et al. (1982) and subsequent 
observations were made by Kawabe (1987) using the Nobeyama Millimeter Array (NMA) 
in 1987. He found that the masers have a doubly-peaked spectrum with a velocity coverage 
of 15 km s _1 and that all the maser spots are concentrated within an area of radius 0."3 
(180 AU). At this position, no free- free radio continuum emission has been detected (Felli 
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et al. 1984; Bally et al. 1983; Hoare et al. 1994). The presence of the H 2 masers, dense 
molecular cloud core and the cold SED indicate that S106 FIR may harbor a protostar. In 
particular, the absence of distinct CO molecular outflows suggests that this protostar is in 
an extremely early state of star formation. 

2. OBSERVATIONS 

Deep integration synthesis observations of the H2O maser line (616 — 523; ^rest = 
22235.077 MHz) toward S106 FIR were performed with the NRAOQ Very Large Array 
(VLA) in its A and BnA configurations on 1996 October 15 and 1997 January 23, 
respectively. The FWHM of the primary beam was 2/0 and the system noise temperatures 
were around 170 K. We observed in the single IF-mode using 128 channels with a total 
bandwidth of 3.125 MHz. The velocity range covered was 42.1 km s _1 with a velocity 
resolution of 0.33 km s _1 at 22.235 GHz (Table |l|). The phase tracking center in the 
A configuration observations was located at the position reported in the previous NMA 
measurements (R.A. 1950 = 20 fe 25 m 32/45, Dec.i 95 o = 37° 12' 50."95; Kawabe 1987). The 
center in the BnA configuration observations was chosen to be at the strongest maser spot 
in the A configuration observations. The integration time per visibility was 8 seconds, and 
the field of view (FOV) was limited to be 37" by a coherency level of 97 % for the longest 
baseline. 

EDITOR: PLACE TABLE ITl HERE. 



x The National Radio Astronomy Observatory (NRAO) is operated by Associated 
Universities, Inc., under cooperative agreement with the National Science Foundation 



The continuum source 2005 + 403 was used as a phase and amplitude calibrator, and 
was observed for 2 minutes during each 15-minute observing cycle. For the first epoch, 
2005 + 403 was also used as a bandpass calibrator by integrating over the run. For the 
second epoch, 3C454.3 was observed as a bandpass calibrator using 3-minute integration at 
the beginning and end of the observations. The flux density scale for the first epoch was 
established using the continuum source 0404 + 768 which was not resolved even in the A 
configuration. The flux density scale for the second epoch was set by the continuum source 
3C286. According to the VLA List of Calibrator Sources, the flux densities of 0404 + 768 
and 3C286 at A = 1.3 cm were 1.0 Jy and 2.5 Jy, respectively, and the accuracy is expected 
to be better than 30%. We made maps of the central regions of 128 mas x 128 mas area in 
the FOV using task IMAGR in the AIPS package with a cell size of 0.5 mas. The resulting 
synthesized beam sizes were 63 mas x 61 mas at P.A. = 85° and 250 mas x 80 mas at P.A. 
= 89° for the first and second epochs, respectively. 

We obtained the absolute positions of the strongest maser spots at the velocity 
channels of V LS r = 6.4 km s' 1 ( R.A.i 950 = 20 fe 25 m 32. s 533 ± 0/00015, 
Dec.i 95 o = 37° 12' 50.950"±0.002") and -13.7 km s" 1 ( R.A. 1950 = 20 ft 25 m 32. s 450±0. s 0006, 
Dec. 1950 = 37° 12' 50.788" ± 0.0025") for the first and second epochs, respectively, before 
a self-calibration procedure. The absolute position accuracies of the strongest spots were 
estimated to be about 0r.a.,dcc ~ 4 mas for the first epoch and or.a. ~ 14 mas, odcc ~ 5 
mas for the second epoch, respectively, taking account of the typical baseline error of the 
array (|Ab| ~ 1 cm) and the angular separation of the phase calibrator from the source 
(| As | = 4.9°). These excellent absolute position accuracies result from the use of the close 
phase calibrator whose coordinates are known with a very high positional accuracy (^ 2 
mas). 
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In order to obtain the relative positions of the remaining maser spots to the strongest 
spots at individual channels and to improve dynamic range of images, we have employed a 
self-calibration procedure using only the phase of the strongest maser spots. The typical 
rms noise levels in the channel maps after the self-calibration were improved to 10.5 mJy 
beam -1 from 28.5 mJy beam -1 and 9.4 mJy beam -1 from 22.5 mJy beam -1 for the first 
and second epochs, respectively. The relative positional accuracies were determined mainly 
by the signal-to-noise ratio (S/N ratio), because the bandpass calibration was done within 
an accuracy of a few degrees across the channels, and because the self-calibration worked 
very well and the residual systematic errors were very small. The relative positional error 
of a point source convolved with a synthesized Gaussian beam is given by 0.45 0/(S/N 



ratio) where 9 is the FWHM of the beam (Reid et al. 1988a; 9 = J9 ma j x 9 min ). For the 
first epoch, we use only the data with S/N ratio higher than 7, resulting in a maximum 
positional error of 5.6 mas (3.4 AU). For the second epoch data, we take the S/N ratio 
higher than 16 which gives the same maximum positional error as that of the first epoch. 
These relative positional errors are comparable to those of the absolute ones except for 
0"r.a. ~ 14 mas of the second epoch. 

EDITOR: PLACE TABLE @ HERE. 

We also monitored the masers irregularly in order to investigate the velocity drifts of 
the maser peaks. The maser emission was observed in 1987 May using the NMA (Kawabe 
1987). Single-dish observations were made with the 45m-telescope of Nobeyama Radio 
Observatory (NRO)Q on 1996 May 5 and 28. Also 5-minute snapshot observations with 



2 Nobeyama Radio Observatory (NRO) is a branch of National Astronomical Observatory, 
an interuniversity research institute operated by the Ministry of Education, Science, and 
Culture of Japan 
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the NRAO, Very Long Baseline Array (VLBA) were conducted on 1996 June 5 as a VLBI 
Space Observatory Program pre-launch survey (Migenes et al. 1998). The observational 
parameters are summarized in Table |2|. 

3. RESULTS 

3.1. Time Variation of the Overall H 2 Maser Spectrum 

Figure [I] shows the H 2 maser spectra obtained by the monitoring observations. The 
intensity scale is shown in Jy unit except for the VLBA snap shot observations because of 
the uncertainties in the flux calibration. The spectra commonly show three major peaks on 
the blueshifted side and a single peak on the redshifted side with respect to the ambient 
cloud velocity (K, ys = —1.1 km s^ 1 ) measured from H 13 CO + J = 1 — observations with the 
NMA (Furuya et al. 1998) and NH 3 (J, K) = (2, 2) observations with the VLA (Mangum 
and Wootten 1994). Distinct maser emission cannot be seen around the ambient cloud 
velocity. The total radial velocity range of the masers is always from Vlsr ~ — 15 to 
12 km s _1 . The blue peaks are commonly seen in the velocity range from Vlsr ~ — 15 to —2 
km s _1 . The red peak is seen from Vlsr ~ +2 to +12 km s _1 , although the peak velocity 
changed from Vlsr = +9 to +6 km s" 1 in the third quarter of 1996. 

EDITOR: PLACE FIGURE HERE. 

We found velocity drifts of the spectral peaks in the spectra of Figure [I]. The velocities 
of the three distinct blue peaks drifted together toward the blue, particularly from May 
1996 till January 1997. The blue peaks at Vlsr = -7.4, - 9.4 and - 12.1 km s" 1 
drifted to Vlsr = —7.8, — 9.8 and — 13.7 km s _1 , respectively, during the two epochs 
of VLA mapping. In addition, we found similar peak velocity drifts on the red side at 
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Vlsr — +9 km s -1 from May 1996 until June 1996 and at Vlsr — +6 km s _1 during the 
VLA observations. The red peak at Vlsr — 9.2 km s _1 drifted to Vlsr = 9.5 km s _1 in 
about one month and the peak at Vlsr = 6.4 km s" 1 drifted to Vlsr = 6.7 km s _1 in about 
three months. Here, we fit Gaussian shapes to the spectra in order to determine the peak 
velocities. These peak velocity drifts can be clearly seen in Figure which shows the peak 
velocities as a function of the observing date. 

The important point which we found from Figure |2| is that individual peaks seem 
to have their own constant accelerations along the line of sight. For the blue peaks, the 
accelerations of the velocity drifts are estimated to be —2.0, — 1.1 and — 1.1 km s~ 1 year^ 1 
for the peaks 1, 2 and 3, respectively; and their mean is —1.4 km s _1 year _1 . For the 
red peaks, the peak at Vlsr ~ 9 km s _1 showed a velocity drift with an acceleration of 
+3.9 km s _1 year _1 on May 1996. The peak at Vlsr ~ 6 km s" 1 showed a velocity drift with 
an acceleration of +1.2 km s~ 1 year~ 1 during the VLA observations. These velocity drifts 
are summarized in Table |3|. 

EDITOR: PLACE FIGURE | HERE. 
EDITOR: PLACE TABLE | HERE. 



3.2. Spatial Structures of Maser Spots 

Figure |3] shows the relative peak positions of the H 2 maser spots in channels of 
0.66 km s _1 resolution, achieved by combining together two adjacent individual channels, 
for the two epochs of our VLA observations. We plotted maser spots with S/N ratios higher 
than 7 and 16 for the first and and second epochs, respectively. Map origins plotted by 
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double-circles show the positions of the phase referencing maser spots. A spectrum of the 
masers integrated over the 128 mas x 128 mas area is also shown with a velocity resolution 
of 0.33 km s _1 at the upper right corner of each panel. Even at the highest resolution of 
0.33 km s _1 , significant differences cannot be seen in the position and position- velocity 
maps when compared with Figures § and ^ at a resolution of 0.66 km s _1 . On the other 
hand, 0.99 km s _1 resolution exceeds the typical line width of the spectral peaks and causes 
a serious degradation on the definition of maser components (see §3.4). Thus, we employ a 
resolution of 0.66 km s _1 in the maps to avoid the confusion caused by many points. 

EDITOR: PLACE FIGURE | HERE. 

We found that two clusters of maser spots are located in the center of the submillimeter 
core around S106 FIR and that no maser spots are found outside the 128 mas x 128 mas 
central area. This result is consistent with the previous NMA results ( Kawabe 1987). The 
apparent separation between the intensity-weighted central positions of the two clusters is 
76.2 mas ( 45.7 AU) at P.A. = 72° and 87.5 mas ( 52.5 AU) at P.A. = 70° for the first 
and second epochs, respectively. Each cluster is elongated along a line connecting the two 
clusters. 

3.3. Velocity Structures of Maser Spots 

Two clusters of the maser spots were found toward the Class source SI 06 FIR. 
The western and eastern clusters of maser spots are ~ 5 — 13 km s _1 blueshifted and 
~ 6 — 10 km s~ x redshifted compared to the ambient cloud velocity, respectively. Hereafter, 
we call them blue and red clusters. The velocity coverages of the blue and red clusters are 
approximately symmetric with respect to the ambient cloud velocity. 



-13- 

In order to investigate internal velocity structures within each cluster, we show 
position-velocity maps along two orthogonal lines in Figure f|. The top and bottom panel 
in Figure f| are the position- velocity maps along and perpendicular to a line connecting the 
two clusters. This line is not only parallel to the elongation direction of the clusters but 
also parallel to proper motions which we will describe in §3.5. We determined the line by 
X 2 -fitting and its position angle was measured to be 72° ± 1° for October 1996 and 70° ± 2° 
for January 1997. We adopt a mean of P. A. = 71°. In Figure [|, filled- and open-circles 
indicate the maser spots for the first and second epochs, respectively, and the vertical 
scales of the position- velocity diagrams are the relative radial velocities with respect to the 
ambient cloud velocity. 

EDITOR: PLACE FIGURE | HERE. 

As is obvious from the top panel of Figure ^, there exist systematic velocity gradients 
within each cluster along the line at P. A. = 71°, particularly in the blue cluster. The 
absolute value of velocity increases with increasing the distance from the middle of the two 
clusters. Mean velocity gradients over the two epochs were —0.74 km s _1 AU _1 for the blue 
cluster and —0.36 km s _1 AU _1 for the red, determined by x 2 -fitting. The velocity gradient 
for the blue cluster is twice as steep as that of the red cluster. The bottom panel of Figure 
H shows the position-velocity diagram along the line at P.A.= —19°. All of the maser spots 
are distributed in a strip of 14 mas (8.4 AU) width centered at the P.A.= 71° line. No clear 
velocity gradient can be seen along the line at P.A.= —19° in each cluster. 

3.4. Identification of Maser Components 

In this subsection, we will identify groups of spatially localized maser spots 
corresponding to the major peaks in the spectra at each epoch. Then we will identify the 
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maser components across the two epochs. 

In order to identify maser components in each epoch, we use the VLA images with 
velocity resolution of 0.66 km s _1 (Figure |3[) and the spectra of the masers (Figures [I] and 
|3|). In the spectra of the masers, we identified three distinct velocity peaks (Vlsr ~ ~ • 7, — 9 
and —12 km s _1 ) on the blue side and a single velocity peak (Vlsr ~ 6.0 km s _1 ) on the 
red side. Here, we call the peaks at Vlsr ~ — 7 km s _1 as Al and Bl on the first and 
second epochs, respectively. Similarly, we label the peaks at Vlsr ~ — 9 and —12 km s _1 
as shown in Tables |] and |5|. In the upper right corner of each panel of Figure |3|, the solid 
bars above the velocity axis show the radial velocity ranges of the three blue peaks and one 
red peak which we will identify as maser components. In Figure |5|, we show maps of the 
relative positions of maser spots by displaying the velocity ranges of the spectral peaks in 
pseudocolors. In the upper right corner of each panel, we show the definition of the velocity 
ranges of the peaks. We found that the maser spots marked by individual colors seem to be 
spatially localized. Hence, the individual spectral peaks can be identified with these spatial 
components. We show the relative positions of the components in Tables f| and |5|. 

EDITOR: PLACE FIGURE | HERE. 

Next, we will show correspondence of the maser components between the two epochs. 
The blue Al, A2 and A3 components on the first epoch are considered to be identical to 
the blue Bl, B2 and B3 components on the second epoch, respectively, from the following 
kinematical arguments. First, each pair shows a similar velocity range over the two epochs. 
Second, the peak velocity drifts over ~250 days in Figure [l] strongly suggest the presence 
of three kinematical components on the blueshifted side (i.e., blue 1, 2 and 3), which 
move with their own constant accelerations (Figure |2[). In addition, the red component 
has a positive acceleration similar to the absolute values of the acceleration of the blue 
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components, although the red component was observed only over the two epochs. Third, 
the blue components show similar relative proper motions with respect to the red one, as 
will be shown in §3.5. 

EDITOR: PLACE TABLE | HERE. 
EDITOR: PLACE TABLE | HERE. 



3.5. Relative Proper Motions of Maser Components 

We found an increase in the separations of about 10 mas over the two epochs between 
the red and three blue components in Figure [5[ This fact means that there exist proper 
motions along the line connecting the two clusters. The absolute positions of the phase 
referenced maser spots in the two epochs were used to bring self-calibrated relative position 
maps onto an absolute position basis. This procedure complicates proper motion estimation. 
Because typical absolute positional errors are cr' RA Doc ~ 7 mas for the first epoch and 
a R.A. ~ 15 mas, a^ ec ~ 8 mas for the second epoch, respectively, and are comparable to the 
proper motions, the proper motions are imprecisely determined. 

Nevertheless, we can derive relative proper motions from comparison of the relative 
position maps whose map origins are at the positions of the red spot at Vlsr = 6.4 km s _1 
for the first epoch and at the red spot at Vlsr = 6.7 km s _1 for the second. This 
superposition is based on our interpretation that the two red spots are physically identical 
as described in §3.4. We show a magnified H2O maser spot map of the blue cluster for the 
two epochs in Figure [5[ Dashed-lines with labels denote the identified maser components. 
We now find the relative proper motions of the components over the two epochs. 
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The distances between the red component and the three blue components with labels 
of 1, 2 and 3 increased from 71.7, 73.2 and 83.7 mas to 84.1, 82.1 and 96.3 mas, over the 
two epochs, respectively. The corresponding proper motions are 12.4 ± 5.4, 8.9 ± 4.3 and 
12.6 ± 3.9 mas for the blue components 1, 2 and 3, respectively. The proper motion derived 
for the component 3 was above the 3<r level, although the other two were at the 2a level. 
The directions of the motions are almost parallel to the line connecting the two clusters. 
Table ^ summarizes the changes of the separations and the relative proper motions. The 
intensity weighted mean of the proper motions of the blue components relative to the red 
one is 10.9 ± 4.2 mas (6.5 ± 2.5AU) per 100 days. This can be converted to a relative 
transverse velocity of 113 ± 44 km s _1 along the plane of the sky. 

EDITOR: PLACE TABLE g HERE. 



4. DISCUSSION 

We will discuss three possible models for the origin of the two maser clusters: (1) 
masers associate with a compact jet-like flow, (2) masers originate near each star of a binary 
system, and (3) masers are generated in tangential parts of a rotating disk. We conclude 
that the jet-like flow model will explain our results well. 

4.1. Compact Jet-like Flow Model 
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4-1.1. Evidence for Compact Outflow and Its Physical Properties 

The elongation of the spatial distributions of the maser spots, the presence of the blue 
and red cluster of maser spots and the relative proper motions of the maser components 
along the elongation strongly suggest that the maser emission is excited in an expanding 
jet-like outflow. We assume that the driving source of the outflow, a protostar, is located 
at the middle of the line connecting the two clusters. This hypothesis is supported by the 
fact that the spectral profiles of the masers have blue- and redshifted symmetric peaks 
with respect to the ambient cloud velocity which is thought to be the same as the systemic 
velocity of the central protostar. 

We will discuss three-dimensional geometry of the jet-like flow based on our results. 
First, we can derive an inclination angle of the outflow out of the plane of the sky from 
the proper motions and radial velocities of the maser components. The inclination angle i 
comes from 

tan , = l^mean-Kysl 

where Vmean is a mean peak velocity over the two epochs and [i is the relative proper motion 
for each blue component. We obtained i = 5.7°, 10.3° and 10.2° for the three blueshifted 
components, respectively (see Table |6|); we adopt an intensity weighted mean of i = 9°. So, 
the outflow axis of S106 FIR is nearly in the plane of the sky. 

Next, we can estimate an outflow velocity range using both the radial velocities and the 
relative proper motions. Since symmetric proper motions for the blue and red components 
are suggested by the symmetric maser spectra with respect to the ambient cloud velocity, 
the outflow velocity is estimated to be one half of the transverse velocity described in 
§3.5 and ranges from 45 to 65 km s _1 , corrected for the inclination angle (see Table |6|). 
This derived velocity range is comparable to that reported by Claussen et al. (1998) who 
measured the proper motions of the H 2 masers in the Class I source IRAS 05413—0104. 
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The escape velocity from the protostar of S106 FIR ranges from 10.3 to 20.6 km s _1 at a 
distance of 25 AU from the protostar, taking the stellar mass to be 1.5 - 6 Mq deduced 
from the observed bolometric luminosity range (24Lq <^ Lboi ^ IOOOLq). Hence, the maser 
spots could not be gravitationally bound to the star: this is consistent with the outflow 
model. 

Furthermore, we estimate the acceleration of the outflow from the radial velocity drifts 
of the maser components; we assume that the motions of the maser components are limited 
to lie along the outflow axis. The resulting acceleration ranges from 4 to 15 km s~ 1 year~ 1 
corrected for the inclination angle. We summarize properties of the outflow in Table |7|. 

EDITOR: PLACE TABLE HERE. 



4-1.2. Formation Mechanism of the Jet-like Flow 

We will discuss two acceleration mechanisms for the jet-like flow. The masers could 
be accelerated by the stellar wind from the assumed protostar, or by stellar radiation. We 
conclude that the former mechanism seems most plausible, because the radiation pressure 
appears insufficient to accelerate the jet-like flow up to the acceleration range estimated in 
§4.1.1. 

First, we will discuss the mechanism whereby the masers are accelerated by the stellar 
wind which is commonly thought to accelerate the water masers (e.g., Genzel et al. 1981). 
Here, we consider each identified maser component as a spherical maser cloudlet with a 
radius of a ~ 5 mas (3AU), which is a typical size of the identified components described 
in §3.4, for simplicity. Since the acceleration, aw, is caused by a pressure gradient over 
the cloudlet, we have aw = - • — ~ - • — , where p is a mass density of the cloudlet 
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and Pw is a pressure from the stellar wind. The mass density of the maser cloudlet, p is 
^H 2 m H 2 > where % 2 is the number density of the cloudlet and mn 2 is the weight of molecular 
hydrogen. The wind pressure, P\y at a distance r from the protostar is given by M ^S , 
where M w and v w are the mass loss rate (typically ~ 10 -7 MQyear~ 1 ) and wind velocity 
(typically ~ 200 km s" 1 ) of protostellar jets, respectively (e.g., Edwards, Ray and Mundt 
1994). Thus, we obtain the acceleration from stellar winds as follows 



«w 



1 1 M w v w 



^H 2 m H 2 a 47rr 2 



2 /Q ATT\ /m8 ^^-3 S 



7 5 10ns- 1 r-M Mw U ^ w W 2 5 AU ^ /3 AUW lO* cnr , 

yr VlO- 7 M yr- 1 ; V200 kms" 1 / V r / V a )\ n H2 J 

(2) 

This value seems to be roughly consistent with the acceleration range of 4 ~ 15 km s _1 year -1 
estimated in §4.1.1. Here, we assumed a typical molecular hydrogen density of 10 8 cm" 3 , 
necessary to produce maser emission. If we consider smaller cloudlets, down to maser spot 
sizes of 1 AU (e.g., Reid and Moran 1988b), «w becomes larger. However this does not 
have a serious influence on our discussion. Hence, it is probable that the jet-like flow traced 
by the masers was mainly accelerated by the stellar wind pressure. 

Next, we will discuss the acceleration, a:^, due to stellar radiation pressure. The 
acceleration from radiation pressure is given by K ^°\ ■ If we take the bolometric luminosity, 
Z-boi to range from 24 to IOOOLq and the mass opacity coefficient, k to be 230 cm 2 g _1 
(a typical value in interstellar medium), then the acceleration from the radiation pressure 
ranges from 0.13 to 5.2 km s _1 year" 1 at r = 25 AU from the star. Here, we assume that the 
maser cloudlets are optically thin. If we take an optically thick maser cloudlet, cxr becomes 
smaller compared with the optically thin case. This estimated range is smaller than the 
inclination-corrected acceleration of ~ 4 — 15 km s" 1 year" 1 . Hence, it seems unlikely that 
stellar radiation pressure accelerates the jet-like flow. 
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4.2. Other Models for the Origin of the Masers 

In this subsection, we discuss the remaining two models for the origin of the masers. 
First, we will discuss a model in which the masers originate near each star of a protobinary 
system. Recent infrared and optical observations of pre-main sequence (PMS) stars revealed 
that the frequency of binary systems in PMS stars is comparable to that in main sequence 
stars (^ 50%; e.g., Mathieu 1994: Ghez et al. 1997: Kohler and Leinert 1998: Brandner 
and Kohler 1998). This fact suggests that binarity must be common among protostars. 

In the binary model, the separation of the two maser clusters can be taken as the 
projected binary separation (2a' = 48 AU). This separation is near the median separation 
of binaries (2a' =30 AU from Duquennoy and Mayer 1991). The orbital plane of the binary 
system is suggested to be perpendicular to the plane of the sky because the proper motions 
were parallel to the line connecting the two clusters. This is a suitable configuration for 
the maser amplification. We estimate a rotational velocity of the binary stars in the orbit 



as H.i n = a/(Vlsr — V sys ) 2 + (|/i) 2 ~ 60 km s^ 1 . A binary phase measured from the zero 
projected separation can be estimated from tan0 = |Vlsr — Kys !/(§/•*) for the individual 
maser components. We adopt a mean value of = +10°. 

The binary model has a serious drawback that the system is not in a gravitationally 
bound state. For the binary system to be gravitationally bound, the following large enclosed 
mass is required. 

*— <r < W ~ 56OM ( T £^) (^^f (3) 

Here, a is the half of the binary separation given by a = a'/sin0 ss 140 AU. This enclosed 
mass is equivalent to a reduced mass of the system. However, the mass is considerably 
larger than the reduced mass of the binary system expected from the bolometric luminosity 
range. Therefore, we reject the protobinary model. 
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Next, we discuss a model of a protostellar disk: a single star is located between the 
clusters and the masers are excited in tangential parts of a rotating disk around the star. 
Elmegreen and Morris (1979) originally studied such a model in an ideal case: if there exists 
enough column density of velocity coherent gas along the line of sight in an edge-on rotating 
disk, then H 2 masers can be emitted in tangential parts of the disk. The resulting maser 
spectrum is expected to have a symmetric double-peaked shape with respect to the systemic 
velocity. The terminal velocity of the spectral peaks corresponds to the rotational velocities 
of the maser emitting regions of the disk. Water masers originating in disks were found, for 
the first time, towards the Class I protostar IRAS 00338+6312 with the VLA (Fiebig et 
al. 1996). The maser profile defined a symmetric double-peaked shape with respect to the 
ambient cloud velocity. The velocity structure of the maser spots suggests that the H 2 
masers are associated with an infalling, rotating disk-like envelope whose radius, traced by 
the masers, ranges from 10 to 30 AU. 

The disk model is not plausible in case of S106 FIR, although the H2O maser spectra 
in the source have a symmetric double-peaked shape. In S106 FIR, the disk plane is 
suggested to be nearly edge-on as mentioned in the binary model, and we can estimate an 
enclosed mass of the disk and star in the same way as we did for the binary model. The 
estimated enclosed mass is large, on the order of IOOMq, which cannot be reconciled with 
the bolometric luminosity range. Therefore, we reject the disk model. 

4.3. Future Works to Reveal the Nature of S106 FIR 

The presence of the H 2 masers and the cold SED suggests that the star formation 
has just begun in this source. We stress that no well-developed outflows have been reported 
in S106 FIR (Hayashi et al 1993 and Bachiller 1996). In addition, all of the Class sources 
except S106 FIR are known to possess distinct CO molecular outflows. These facts suggest 
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that S106 FIR is at an evolutionary state just after the onset of outflow activity and the 
formation of well-developed CO molecular outflows usually observed. Do these masers 
represent a very early state in the development of an outflow from a protostar ? 

We have searched for bipolar molecular outflows using the lines of CO J— 1-0 and 
HCO + J— 1-0 with the NMA and the NRO 45m-telescope, respectively. However, no 
large-scale molecular outflow has been found (Furuya et al. 1999). On the other hand, 
H 13 CO+ J=l-0 observations with the NMA and NH 3 (J,K) = (l,l) and (2,2) observations 
with the VLA showed the presence of a dense molecular cloud core around S106 FIR with 
a mass of ~ 1.0 — 1.5Mq and a size of 9000 AU. This dense core is expected to supply 
molecular gas to a large-scale outflow through interaction with a stellar wind or the compact 
jet-like outflow. 

Our work is being continued as further studies of the proper motions of the H2O 
masers with the VLBA. The detailed structure and proper motions of the maser spots will 
be discussed in a forthcoming paper together with a search for large-scale CO outflow with 
the NMA. A search for the exact position of the central star(s) by free-free continuum 
observations will also be attempted. We propose that the central star would be located at 
the middle of the two clusters ( R.A. 1950 = 20 h 25 m 32/49, Dec.1950 = 37° 12' 50."87). 

5. SUMMARY 

High angular resolution VLA observations of the 22 GHz water maser line were 
conducted toward the Class source SI 06 FIR. The main findings of this paper are as 
follows. 

1. We found two clusters of the H 2 maser spots with a separation of 50 AU, located at 
the center of the submillimeter core around the Class protostar S106 FIR. The sizes 
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of the clusters are 50 AUx5 AU. The western and eastern clusters are 8.0 km s _1 
blueshifted and 7.0 km s _1 redshifted, respectively, with respect to the ambient cloud 
velocity. 

2. Three maser components were identified in the blue and a single one in the red 
clusters, respectively. These components are composed of spatially localized maser 
spots within a 10 mas (6 AU) area, and correspond to the major peaks in the overall 
spectrum of the H2O maser line. 

3. We determined both the relative proper motions and peak velocity drifts of the maser 
components. The proper motions are along the line connecting the two clusters and 
are ~ 30 mas year -1 . All the components showed similar peak velocity drifts with 
radial accelerations of ~ 1 km s _1 year _1 . 

4. We conclude from our observations that the masers are associated with a highly 
collimated compact accelerating jet-like flow which could originate from an assumed 
protostar located at the middle of the line connecting the two clusters. 

5. It is possible that SI 06 FIR is an extremely young protostar in an evolutionary 
state just after the onset of outflow activity, because of the absence of distinct CO 
molecular outflows. 
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Table 1: VLA Mapping Observations of H 2 Masers toward S106 FIR 



Epoch 


Config. 


Integ. Time 


Synthesized 
Size 


Beam 
P.A. 


Vel. Coverage 


Vel 


. Resolution 






(hours) 


(mas) 


(cleg) 


(km s _1 ) 




(km s _1 ) 


1996 Oct 15 


A 


3 


63 x 61 


85 


42.1 




0.3 


1997 Jan 23 


BnA 


2 


250 x 80 


89 


42.1 




0.3 



This manuscript was prepared with the AAS IATj^X macros v4.0. 

Fig. 1.— Spectra of the H2O maser emission toward S106 FIR observed with the NMA 
(1987 May 8), the Nobeyama 45m-telescope (1996 May 5 and 28), the VLBA (1996 June 5) 
and the VLA (1996 October 15 and 1997 January 23). The intensity scales are shown in Jy 
unit except for the VLBA observations. The ambient cloud velocity (V^ ys = — 1.1 km s _1 ), 
measured from H 13 CO + J = 1 — and NH3 (J, K) = (2, 2) observations, is also shown. 
Arrows drawn on the spectra show the drifts of the major peaks ( see in text). 
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Fig. 2.— Time variations of the velocities of the major peaks of the H2O maser spectrum 
in S106 FIR on the redshifted (upper panel) and on the blueshifted (lower panel) sides. The 
peak velocities are determined by applying Gaussian profiles to individual peaks. 

Fig. 3.— Peak position maps of the H2O maser spots shown with a velocity resolution of 
0.66 km s" 1 with S/N ratio higher than 7 for (a) 1996 October 15 observations and higher 
than 16 for (b) 1997 January 23. The map origins shown by double-circles are defined by 
the positions of the phase referencing maser spots at the Vlsr = 6.4 km s _1 component on 
1996 October and at the Vlsr — —13.7 km s~ x component on 1997 January. The spectrum 
of the masers integrated over the mapping region is also shown at the upper right corner of 
each panel with a velocity resolution of 0.33 km s _1 . Solid bars above the velocity axis in 
the spectra indicate the ranges of the identified velocity components. 

Fig. 4.— Upper and lower panels: superposed position-velocity maps of the maser spots 
on 1996 October 15 observations (filled circles) and on 1997 January 23 (open circles). The 
upper and lower panels show the position - velocity maps along the line at P.A.= 71° and 
P.A.= —19°, respectively. In the ordinate of the position- velocity maps, the ambient cloud 
velocity (V sys — — 1.1 km s _1 ) is subtracted from the radial velocities of the masers spots. 
Center panel: superposed maps of the maser spots over the two epochs. The two orthogonal 
lines indicate the lines of P.A.= 71° and P.A.= —19°, along which we made the position- 
velocity maps. 

Fig. 5.— Peak position maps of identified maser components shown by pseudocolor. The 
definitions of the components are indicated in the upper right corner of each panel. The 
map origins are the position of the 6.4 km s _1 spot for the first epoch, and the position of 
the 6.7 km s _1 spot for the second epoch. The phase referencing maser spots in the maps 
are shown by white-circles. 
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Fig. 6.— Superposed map of the maser spots at the two epochs for the blue cluster. The 
positional accuracies are the same as in Figure 3, and the map origins are the same as in 
Figure 5. Filled and open circles represent 1996 October and 1997 January observations, 
respectively. Dashed lines with labels indicate the identified maser components. 



Table 2: Summary of H 2 Maser Observations toward S106 FIR 



Beam Size 



Date Telescope Primary Synthesized Vel. Resolution Sensitivity 3 

(") (mas) (km s _1 ) (mJy) 



1987 May 8 


NMA 


330 


5000 


0.26 


~100 


1996 May 5 


NRO 45m 


74 




0.50 


180 


1996 May 28 


NRO 45m 


74 




0.50 


240 


1996 Jun 5 


VLBA b 


120 


« 30 c 


0.21 


220 


1996 Oct 15 


VLA-A 


120 


63 x 61 


0.33 


10 


1997 Jan 23 


VLA-BnA 


120 


250 x 79 


0.33 


10 















"Rms noise level per velocity channel. 

6 Performed as a pre-launch survey of VLBI Space Observatory Program (VSOP). 

c The size is estimated from the fringe spacing of the shortest baseline, between Los Alamos 
and Pie Town in the VLBA. 
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Table 3: Radial Velocity Drifts of H 2 Q Maser Peaks in S106 FIR 



Maser Component 


Radial Velocity 


Radial Acceleration 


15 Oct 1996 


23 Jan 1997 




(km s _1 ) 


(km s _1 ) 


(km 


s l yr v ) 


Red 


6.4 


6.7 




+ 1.2 


Blue 1 


-7.4 


-7.8 




-2.0 


Blue 2 


-9.4 


-9.8 




-1.1 


Blue 3 


-12.1 


-13.7 




-1.1 


Blue (mean) 








-1.4 a 













"Intensity weighted mean 



-32- 



Table 4: Relative Positions of H 2 Maser Components in S106 FIR on the First 
Epoch (1996 Oct. 15) 



Maser Components 



Position Offset 3 



Total Flux Density 



Peak Velocity 

(km s _1 ) 



Name 



A R.A. 



mas 



A Dec. 



mas 



(Jy) 



6.4 b 




0.0 


0.0 


1.9 


-7.4 


Al 


68.9 ±1.6 


-19.7 ±1.6 


0.3 


-9.4 


A2 


70.5 ±0.4 


-19.7 ±0.4 


1.1 


-12.1 


A3 


81.3 ±0.4 


-19.8 ±0.4 


1.3 













"Intensity weighted mean position 
6 Map origin 
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Table 5: Relative Positions of H 2 Maser Components in S106 FIR on the Second 
Epoch (1997 Jan. 23) 



Maser Components 



Position Offset 3 



Total Flux Density 



Peak Velocity 
(km s _1 ) 



Name 



A R.A. 



mas 



A Dec. 



mas 



(Jy) 



6.7 




-93.8 ±2.8 


22.0 ±2.7 


1.1 


-7.8 


Bi 


-12.3 ± 1.6 


1.09 ±1.0 


1.6 


-9.8 


B2 


-1.7 ±3.4 


-13.8 ±2.6 


2.2 


-13.7 b 


B3 


0.0 


0.0 


4.0 













"Intensity weighted mean position 
6 Map origin 
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Table 6: Relative Proper Motions of H 2 Maser Components in S106 FIR 



Name Distance from Red Comp. /i a V mcan i h V 3 v c 

15 Oct 1996 23 Jan 1997 

(mas) (mas) (mas) (km s _1 ) (degree) (km s _1 ) 



Blue 1 (A1->B1) 
Blue 2 (A2^B2) 
Blue 3 (A3^B3) 



71.7±2.3 


84.1±4.9 


12.4±5.4 


-7.6 


c 7+3.6 
°- '-1.9 


64±J:g 


73.2±0.6 


82.1±4.3 


8.9±4.3 


-9.6 


10.3^ 6 2 


46±iS 


83.7±0.6 


96.3±3.9 


12.6±3.9 


-12.9 


IQ.2^3 5 


65±S:i 



"Relative proper motions measured from the red component. 
^Inclination angle defined by arctan ( ' mcan ~ ml \ where V sys — ~~ 11 km s" 1 . 
Three-dimensional velocity obtained assuming symmetric proper motions for the blue and 
red components. 
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Table 7: Physical Properties of Compact Jet-like Flow in S106 FIR 

Size (AU) 50x5 

P. A. (degrees) 70 

Inclination (degrees) 9 

Outflow Velocity 21 (km s" 1 ) 45~65 

Outflow Acceleration 13 (km s _1 year -1 ) 4~15 

"Traced by the masers. 

6 Acceleration, traced by the masers, corrected for the inclination angle of 9° assuming that 
the motions of the maser components are along the outflow axis. 
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